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Motivation 
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Schematic figure: 


Temperature 7 


Nuclear Superfluid us 


(Figure is adopted from Fukushima, Hatsuda, Rept. Prog. Phys. 74 (2011) 014001) 
@ Finite temperatures --- Lattice simulation works well ! 


ə Finite density --- No ab initio calculations. 
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Sign problem of finite-density QCD 


QCD 
Fundamental theory for quarks and gluons 


Neutron star 


ə Cold and dense nuclear matter 


@ 2Msun neutron star (2010) 


Neutron star merger 
ə Gravitational-wave observations (2016~) (image from NASA) 


Path integral for finite-density QCD; 


Zoco (T, 1) = f DA Det(P(A, u) +m) exp Sul), 


quark gluon 


Sign problem: Det(D(A, uy) +m) #0 at uq # 0. 
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Nonperturbative methods in QFT 


Numerical LATTICE QFT 
ə Manifest Gauge invariance, Broken spacetime symmetries 
ə Need Monte Carlo integral. Suffer from the sign problem 


Functional method 
ə Do not perform the path integral explicitly 
ə Derive the consistency functional relation and Solve it instead. 


ə Manifest spacetime symmetry 
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Example: Dyson-Schwinger equation 
Quantum version of equation of motion; 


(SH) (a1) --- O(an)) = 0. 


Diagrammatically, 


Other examples 
oe 2PI method 


ə Functional Renormalization Group 


Question 
Do not these functional methods suffer from the sign problem? | 
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Problem 
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Basic procedures of Functional methods 
BCS theory of superfluidity: 


S =4* (3; — V? — u) + Ary + diet) + A?/ (29) 


1. Compute free energy F(A): Mean-field approximation 


A? d’k iwn + k? — u —A 
FA) = ae, (m; In det ( IA ie ik, ) 


2. Solve the equation of motion 
F'(A) = 0 is the BCS gap equation. 


3. Pick up the solution A, that minimizes F 
e For T > T,, A, = 0 is the unique real solution. 
e For T < Te, U(1) symmetry is spontaneously broken, A, # 0. 
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Effective potential and the sign problem 


Path-integral definition of the constrained free energy is 
exp (-BVF(A)) = | Doexp (-SIo)) (ALE - A) 


If the microscopic action S[@] is real, then F(A) is real. 
=> All the basic procedures can be justified by thermodynamics. 


For complex S|], F(A) becomes complex. 


Problem 
e What is the physical meaning of complex F(A)? 
ə Which space must be considered for A? 


ə If there are several solutions of F'(A) = 0, which one should be 
selected? 
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Solution 


Solution 
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Meaning of complex F(A) 


Problem 
e What is the physical meaning of complex F(A)? | 


It is impossible to give the physical meaning of F(A) at each value 
of the background field A € R. 


Thermodynamic free energy Finermo. 
exp (BV Finermo.) = | Pe exp (—S|d]) 5 
can be expressed using F'(A) as 


exp (AV Fihermo.) = [aa exp (-BV F(A)) z 


This equation is the physical meaning of F(A). 
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Appropriate space and selection of A 


Problem 
@ Which space must be considered for A? 


ə If there are several solutions of F'(A) = 0, which one should be 
selected? 


Since we are interested in thermodynamics, evaluate 


exp OV Finermo.) = [aa exp (-BV F(A)) : 


in the limit V — oo. 


Use the steepest descent method (Lefschetz-thimble method) for this 
purpose! 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = ri (= + ar) 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a= 1) 
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Rewrite the Airy integral with Lefschetz thimbles 
There exists two Lefschetz thimbles J, (a = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No f. s pi (= +a) 


Ng: intersection number of the steepest ascent contour Ke and R. 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles 7,: (classical eom F’(z,) = 0) 


i dx ee) = > (Ke; R) | d”ze FD, 


[02 


Js are called Lefschetz thimbles, and Im[S] is constant on it: 


in se a dzi(t) _ eo) 


t——-00 dt Ozi 


I= 420) 


(Ko, R): intersection numbers of duals K, and R” 


(Ko = {2(0)|2(00) = 2, }). 
[Witten, arXiv:1001.2933, 1009.6032 in semiclassical analysis of complex CS theory] 
[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc. in sign problem of LFT] 
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Constrained free energy and Lefschetz thimbles 


Let's go back to our problem, then 
exp (BV Finermo.) = f dA exp (-BV F(A)) 
R 


= ERK) f dA exp(-BVF(A)). 


o 


One can now estimate RHS at saddle points, and as V — oo 


exp(—8V Finermo) ~ X (R, Ko) exp(—BVF(A,)). 


oO 
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Procedures for functional method with sign problem 


1. Compute the constrained free energy F'(A) 
This part is the same. 
2. Solve the equation of motion “for complex A” 


Solve F’(A) = 0 in complexified space, and denote the complex 
solutions as {A,},. 


2’. Solve the gradient flow to construct J, and Ko 


dt \dA 
3. Select A, that minimizes Re(F(A)) with (R, Ko} # 0. 
Since F = F(A,), F(A) must be real for consistency. 


dA OF 
This step is completely new. One has to solve — = ( ) 
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Application 
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Application 


Polyakov-loop effective model 


What can be a model of low-energy QCD? 
e Confinement (Center symmetry): Polyakov loop 


1 B 
{3 = -tr P exp (i ara.) 
3 0 


ə Mass generation (Chiral symmetry): chiral condensate 
(vd) 


Ginzburg-Landau theory must contain (3), (£3), and (yy). 
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Sign problem of effective model at finite densities 


Quark determinant at finite u with nontrivial Polyakov loop is 
complex. (Dumitru, Pisarski, Zschiesche, PRD72(2005)065008, 


Fukushima, Hidaka, PRD75 (2007) 036002) 


Effective potential for heavy-dense quark (simplified version) is 
3—1 
Ver = „ED GRACO 02) + e tai, 02)) ’ 
where 


(eer?) + el(— 61462) + e 72) 


en = (e i(91+82) | e i( 91402) Leif). 


At u #0, Veg takes complex values. 
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Application 


Charge conjugation 


Free energy is given by 


2 01 +302 . > 01 — 305 
sin 5) 


exp(-F) = fanan, sin? 0, sin 
x exp(— Ver (01, 02)). 
Charge conjugation interchanges (3 and (z, which gives 
Ver (01, 02) = Ve (61, —02) 


Since the charge conjugate pair is summed up, exp(—F’) is real. 
(Fukushima, Hidaka, PRD75 (2007) 036002, YT, Nishimura, Kashiwa, PRD91 (2015) 101701) 
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Application 


CK transformation and Gradient flow 


Assume that F(A) satisfies, under “charge conjugation” C 
(C=C =C"), 


F(A)=F(C-A) for AER". 


Define the CK transformation for A(€ C) as 
CK(A)=C-A. 


(Nishimura, Ogilvie, Pangeni, PRD90 (2014) 045039, PRD91 (2015) 054004) 
This is compatible with the gradient flow, i.e. denote A(t) as a 


solution of 
dA(t) = OF (A) 
dt OA i 


then CK(A(t)) satisfies the same equation. 
(YT, Nishimura, Kashiwa, PRD91 (2015) 101701) 
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Application 


Refinement of Lefschetz-thimble formula with CK 
Decomposing the set of complex saddle points F’(A) = 0 as 


Bp = {Ay =CK(A,)}, Ea = {Ac |Im(F(A,)) 2 0}. 


+ (RK) [ ange d”A exp(—F(A)). 


TEX 


In this formula, each term is real because each integral sums CK 
conjugate pairs. 
(YT, Nishimura, Kashiwa, PRD91 (2015) 101701) 
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Application 


Gradient flow of Polyakov-loop model 


ee ee ee ee eT 


X 
` 
` 
` 
` 
` 
` 
` 
` 


o 2 


: Re(z1) 
Figure: Flow at h = 0.1 and u = 2 in the Re(z1)-Im(z2) plane. 

The black blob: a saddle point. 

The red solid line: Lefschetz thimble 7. 

The green dashed line: its dual K. 


(YT, Nishimura, Kashiwa, PRD91 (2015) 101701) 
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Application 


Saddle point approximation at finite density 


The saddle point approximation can now be performed at 


CK-invariant saddle. 


Polyakov-loop phases (21, z2) takes complex values, so that 


(£), (0) € 


Since Im(z2) < 0 and 


L~ 


w| = 


we can confirm that 
L>e. 


(YT, Nishimura, Kashiwa, PRD91 (2015) 101701) 


R. 


(2e'” cos 6, + e 7%), 
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Conclusion 


Conclusion 
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Conclusion 


Conclusion 


ə Functional methods also suffer from the sign problem (but in a 
different way from Numerical Lattice Field Theory). 


ə Using Lefschetz thimbles, we can derive a systematic procedure 
for this situation. 

ə Extension of charge conjugation ensures the computed free 
energy becomes real. 


ə Polyakov-loop effective model is a good object to be 
reconsidered in this framework. 
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